Local structural variation as source of magnetic moment reduction in BaFe 2 As 2 
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We report time-of-flight neutron powder diffraction results on stoichiometric BaFe2As2. Pair 
distribution function analysis shows that the orthorhombic distortion in the a — b plane at short 
distances are significantly different from the average lattice distortion, indicating local variations in 
the lattice at the short-range. We propose that this local variation reflects a high density of nano- 
twins, short-ranged structures which locally affect the magnetic alignment. This results suggests 
that the discrepancies between the observed and calculated magnetic moments in BaFe2As2 arise 
partly from short-ranged variation of the lattice in the a — b plane. 



I. INTRODUCTION 

Since the discovery four years ago of high-temperature 
superconductivity in iron based compounds (FeSC)i, 
much progress has been made to identify the relevant 
physics of these materials. Nevertheless, many questions 
about the superconducting and magnetic properties re- 
main, many of these dealing directly with the questions 
of the influence of structural transitions to the electronic 
properties of the system. We have used pair distribu- 
tion function analysis to describe the local atomic struc- 
ture of BaFe2As2, a parent compound to several fami- 
lies of FeSC. In this Letter, we report the presence of 
local atomic disorder in BaFe2As2, and discuss the im- 
pact of the disorder on the magnetic properties. By com- 
paring magnetic properties measured with spectroscopic 
probes, calculations, and neutron diffraction, we under- 
score the importance of the structural contribution to the 
electronic properties of the system. 

BaFe2As2 is a tetragonal paramagnet above 140 K 
at ambient pressure. Below T m 140 K, the symme- 
try of BaFe2As2 is lowered to an antiferromagnetic (AF) 
orthorhombic state [Fig. [TJi]. This distortion changes 
the planar layer of iron atoms from squares to rectan- 
gles by expanding along the a-axis and contracting along 
the 6-axis by 0.01 A- 2 -^. Superconductivity is induced in 
BaFe2As2 through suppression of the magnetic ground 
state via a number of mechanism (See Refi^ for review) , 
and the magnetic and superconducting behaviors of FeSC 
have clear dependence on the structural properties 4 - - — . In 
the local moment picture, this small orthorhombic distor- 
tion induces large anisotropy in the magnetic exchange 
constants as determined from spin wave fittings, with 
AF Ji a and FM Ju, having an opposite sign [Fig. [TJd] 7 ' 8 . 
Pair distribution function (PDF) analysis of BaFe2As2 
and related compounds by x-rays shows no evidence 
local distortion in the K-doped regime^, but neutron 
analysis shows the presence of spin-ladder structures in 
chalcogenidesi^, and an eventual ambiguity in the struc- 
tural assignment of the Fe(Se,Te) class of materials^. 



Density functional theory (DFT) calculations pre- 
dict values for the magnetic moment from 1.75 to 
2.4 /ig/Fe^ ' 12 ' 13 and show that the local Hund's coupling 
is primarily responsible for large magnetic moment. 8 The 
value of the local moment observed from core electron 
spectroscopy is 2.1 (ib for SrFe2As2^ On the other 
hand, the magnetic moment of BaFe2As2 determined 
from neutron diffraction measurements range from 0.87 
to 1.04 /ig/Fe ; 15 ' 16 and sum-rule integrations over the 
spin excitation spectrum yield 1.04 fis- [13] A*SR and 
Mossbauer studies of BaFc2 As2 show an ordered iron mo- 
ment of Q.4/iB.@and[I|| 

Part of the discrepancy between the core electron 
spectroscopy measurement and other spectroscopic tech- 
niques of nearly a factor of 2 could be dynamic, because 
the time scale of spectroscopic measurement is very short, 
and the DFT calculations do not include quantum fluc- 
tuations. However, it is possible that the discrepancy 
originates also from the fact that the spectroscopic mea- 
surement measures local moment, whereas the diffraction 
measures the average moment. This is the point we in- 
vestigate in this work. 

We have conducted a study of the BaFe2As2 system 
using neutron powder diffraction and pair distribution 
function (PDF) analysis. Information about local distor- 
tion and local atomic correlation is present in the high Q 
diffraction data, for instance as diffuse scattering. The 
PDF is the Fourier transform of S(Q), and retains all of 
the information regarding local structural distortions and 
atomic correlations, allowing one to incorporate this in- 
formation into a comprehensive model of the local struc- 
ture of a syste m 19 i 20 . 



II. EXPERIMENT 

Neutron timc-of-flight powder diffraction experiments 
were conducted at the Manuel Lujan Neutron Scat- 
tering Center at Los Alamos National Laboratory us- 
ing the NPDF powder diffractometer—. A sample of 
BaFe2As2 was prepared using powder synthesis methods 
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FIG. 1. (color online) (a) Orthorhombic crystal structure and 
spin alignment of BaFe2As2 as suggested by neutron scatter- 
ing experiments, (b) Magnitude of experimentally determined 
exchange constants decreases with the relationship Ji a > Ji 
> Jib- (c) Rietveld refinement pattern for BaFe2As2 at 30K. 
(d) PDF data compared to the Rietveld structural model. 
The PDF for the Rietveld obtained structural model deviates 
most substantially at the shortest range, exhibiting a clear 
crossover at fO A. Refinement parameters are summarized in 
Table U 



FIG. 2. color online) Comparison of short-range (2-4.2 A) 
PDF data to the (a) Rietveld and (b) refined PDF struc- 
tural models with interatomic distances labeled. The Rietveld 
structural model deviates in peak positions at the Fe-As, Fe- 
Fe, and As- As distances, while the peak heights indicating the 
coordination in the local environment are not accurate. In b) 
the lattice parameters and isotropic thermal parameters are 
free parameters to the fit, improving the fit markedly. The * 
in (b) marks the distance of the tetragonal a lattice parame- 
ter. 



(Tat = 138 K). Samples were sealed in a vanadium sam- 
ple can in a helium glovebox, and placed into a cryostat 
for temperature controlled studies. The diffraction data 
were analyzed by Rietveld refinement as implemented in 
GSAS^, and the results from the Rietveld refinement 
are shown in Fig. [Ib. Refinement parameters for Ri- 
etveld analysis included lattice parameters, absorption, 
isotropic atomic displacement parameters (ADP), and 
the z position for the As atoms. The background was 
modeled using a Chebyshev polynomial. 

PDF data were produced using standard methods^ 
and quantitative analysis was done with PDFgui 24 . 
Fourier transforms were terminated at 35 A -1 . The Ri- 
etveld structural model was the starting model for the 
PDF analysis. The measured PDF was found to deviate 
substantially from the PDF calculated for the model ob- 
tained by the Rietveld refinement [Fig. [T}i]. Agreement 
is poorest for the short-range (2-4.2 A) [Fig. UK] and im- 
proved at long-range (10-20 A) [Fig. [Tp!]. Allowing lattice 
parameters and isotropic ADPs to be free parameters to 
the fit, we find that the fit in the short-range is markedly 
improved [Fig. ^jp] and that different lattice parameter 
values describe the system for the different fitting ranges. 
Consequently the Fe-Fe distance changes as a function 
of fitting range [Fig. [3]. At the same time the isotropic 
ADPs of Fe and As [Table Q] exhibit a monotonic increase 
in value as a function of fitting range. 



We characterize the behavior of the Fe-Fe distance in 
the a-b plane, i?Fc-Fo, a and -Rf c _f c , b through a quantity 
orthorhombicity, defined as: O = Rfo F ° r a " flF °- Fc ' b . The 
value of the O is maximal at the shortest fitting range at 
1.38%. For comparison, the value determined from the 
Rietveld refinement is 0.78%, representing a 60% increase 
in O at the short-range. 



III. DISCUSSION 

The poor fit of the Rietveld derived PDF and resul- 
tant increase in O at the short-range shows that the lo- 
cal structure of BaFe2As2 seen by PDF is different from 
the average structure determined by the Rietveld anal- 
ysis. The structural heterogeneity is also seen in the 
monotonically increasing isotropic atomic displacement 
factors. Apparently the iron plaquette is locally more dis- 
torted than the average structure [Fig. [3]). Conversely, 
the short-range structure which fits the PDF at short 
distances does not describe the medium- or long-range 
structure [Fig. [4] , showing that the local structure is dis- 
tinct from both the medium- and long-range structure. 

We note that the previous X-ray study (Ref. [^) does 
not report local distortions in BaFe2As2, preferring to 
treat the system with the inclusion of anisotropic ther- 
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TABLE I. PDF and Rietveld refinement results for BaFe2As2 at 30K; u values are isotropic atomic displacement parameters, 
calculated from the trace of the atomic displacement tensor. 



Range (A) 



a (A) 



b(A) 



;(A) 



u-Ba (A 2 



u-Fe (A 2 



u-As 



'-) 



2.0-4.2 
4.0-10.0 
10.0-20.0 
Rietveld 



5.63(2) 
5.620(5) 
5.617(1) 
5.6065(2) 



5.56(2) 
5.566(5) 
5.570(1) 
5.5627(7) 



12.925(4) 
12.953(6) 
12.951(3) 
12.9284(4) 



0.0026(3) 
0.0017(6) 
0.0032(6) 
0.0025(5) 



0.0017(3) 
0.0022(2) 
0.0034(2) 
0.00343(29) 



0.0008(3) 
0.0015(3) 
0.0016(3) 
0.00330(32) 



2.84 



2.78- 



2.76- 
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FIG. 3. (color online) Iron bond lengths from fit PDF for 
BaFe2As2 at 30K. The length of the long side of the plaquette 
gradually decreases as a function of fitting range reflecting the 
degree of short-range distortion in BaFe2As2. 



mal factors. We believe the discrepancies between these 
reports arises from the difference in experimental condi- 
tions: their Q range is limited to 25 whereas our 
neutron measurements reach a higher Q range, up to 35 
A" 1 providing better real-space resolution. 

We explain the difference between the local and global 
Fe-Fe distances using a nanoscale twinning model. In 
this model the local structure has orthorhombic distor- 
tion with O = 1.38%. If the volumes of the nano- twins 
in two orientations are equal, the macroscopic distortion 
will be zero. In our model the volumes of the twins are 
not equal, resulting in a reduced, but non-zero, strain of 
O = 0.78%. An example of such nano- twins is shown 
in Fig. [5] Regular structural modulation as shown in 
Fig. [5] will produce superlattice diffraction. Absence of 
superlattice diffraction with BaFe2As2 suggests that the 
modulation is not regular, but is random, and we note 
that random modulation will not result in broadening of 
diffraction peaks 2 ^. 

The presence of such nano-twins should have a pro- 
found effect on the magnetic structure, because the mag- 
netic exchange parameters are positive (AF) along the 
direction of longer Fe-Fe distance (a axis) , and negative 
(FM) along the direction of shorter Fe-Fe separation (b 
axis)2& Therefore the direction of the AF spin align- 
ment should largely follow the orientation of the twins. 
The simplified model shown in Fig. [5] shows that the 



realignment of the spins along the nano-twin boundaries 
results in nano-scale AF domains and weakening of the 
magnetic diffraction peaks observed by neutron diffrac- 
tion. If, as is shown in the figure, the nano-twin with the 
volume fraction of 20% reduces the average orthorhom- 
bic distortion by 40%, the AF sublattice moment should 
also decrease by 40%. Again, the model shown in Fig. [5] 
is overly simplified and the AF domain walls may have 
some width, resulting in smoother local variations in the 
spin orientation. However, it is obvious that the com- 
plex behaviors of the nano-twins could easily result in a 
significant moment reduction. 

We note also that the value of the moment obtained 
from Mossbauer— is consistent with the value obtained by 
/iSR 18 , a value of 0.4 ^b, which introduces a discrepancy 
of a factor 2 between the slower timescale methods and 
the value of the moment as measured by neutron diffrac- 
tion. We can infer from the reported information that 
the expected error on these values would still place the 
value of the moment less than that reported by neutron 
diffraction. Further, /xSR studies report the presence of 
inequivalent magnetic sites in the unit cell of BaFe2As2 18 , 
raising a question of the true distribution of the spin 
density in the system, and emphasizing the important 
dynamic aspects of the magnetism in this system. 

The nano-twin model is consistent with results in sev- 
eral published works. Above Tjv, the system exhibits 
broad quasi-elastic signals at the AF wave vector in in- 
elastic neutron studies^^, indicating retention of strong 
in-plane medium-range spin order in the paramagnetic 
state in spite of a loss of inter-planar correlation. As 
the temperature of the system is lowered toward Tjv, the 
spatial range of spin correlation increases and the fluc- 
tuation frequency is lowered. Thus spin correlation cou- 
ples with the lattice, resulting in local distortion. This 
effect is represented by the softening of the transverse 
phonon mode near T N ZL The softening occurs over a 
finite range of inelastic momentum transfers q, so the 
correlation length of the fluctuations never diverges as in 
the standard second-order phase transition. Indeed this 
transition has first-order signatur o 26 i 28 ' 29 , suggesting the 
local nature of the transition. Further, the resulting in- 
crease in the magneto-structural transition temperature 
by compressing the orthorhombic 6-axis^ is naturally ex- 
plained by this model. Applying uniaxial strain to the 
6-axis induces partial de-twinning, resulting in the in- 
crease in the volume fraction of the twin with one sense 
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FIG. 4. (color online) Application of the short-range struc- 
ture applied to the medium- (4-10 A) range structure. The 
short-range structural model fails to describe the medium- 
range structure, indicating that the short-range structure of 
BaFe2As2 is distinct. 



and the intensity of the magnetic diffraction peak. 

Additionally, systematic differences in the mag- 
netic moment from neutron diffraction are observed in 
BaFe2As2 as a function of cooling temperature. As ob- 
served in Ref. [U the level of hysteresis in the mag- 
netic transition is highly sample dependent, in addi- 
tion to depending on the rate of cooling used in the 
experimen t 15 ' 16 ' 28 , highlighting the instability of the for- 
mation of the magnetic structure in these systems. 



IV. CONCLUSION 

To conclude, we present evidence for a local structural 
variation of the BaFe2As2 system, in the form of nano- 
twins. The presence of the nano-twins introduces nano- 
scale AF domains in magnetic order in the BaFe2As2 
system, lowering the average static magnetic moment of 
the system seen by neutron diffraction. Observations are 
made in the as-grown compound BaFe2As2, absent any 
distortion that might be expected from chemical substi- 
tution or alloying, showing that the system exhibits a 
fundamental structural and magnetic frustration. By fo- 
cusing on the parent compound we decoupled the effect 
of alloying and substitutional disorder from the interpre- 
tation. But there is ubiquitous evidence of similar lo- 
cal distortions across many FeSC compounds and this 
will be discussed elsewhere. This result addresses an 



outstanding question in the study of high-temperature 
iron-based superconductors, explaining in part that the 
source of the discrepancy between moments measured by 
the spectroscopic probes, DFT calculations, and neutron 
diffraction is intrinsic structural distortion in the pnictide 
compounds, 
b 




FIG. 5. (color online) Suggested nano-twins of BaFe2As2 with 
corresponding alteration of magnetic alignment. Nano-twins 
locally reorient the iron plaquette, forcing the spins to realign 
into nano-scale AF domains. This results in an overall reduc- 
tion of the magnetic moment of the system as measured by 
neutron diffraction. 
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